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Abstract 

We present a numerical study of the properties of the Fixed Point 
lattice Dirac operator in the Schwinger model. We verify the theoretical 
bounds on the spectrum, the existence of exact zero modes with definite 
chirality, and the Index Theorem. We show by explicit computation that 
it is possible to find an accurate approximation to the Fixed Point Dirac 
operator containing only very local couplings. 
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1 Introduction 



Lattice QCD is a tool to perform (in principle exact) non-perturbative com- 
putations of physical quantities in problems where the strong interaction effects 
arc important, by means of numerical simulations. Yet, available computer re- 
sources do not allow simulations on lattices extended over a large enough phys- 
ical size (to keep the finite size effects under control) and with small enough 
lattice spacing (to make the non-physical cut off dependence negligible). 

In the pure gauge sector, together with the contamination of the physical 
results by lattice artifacts, we face up to an unpleasant situation: the difficulties 
in studying the topological effects, a genuinely non-perturbative problem. At 
the classical level, it is possible to define a lattice topological charge operator 
which, as in the continuum, assumes integer values on all lattice configurations, 
the Liischer's jl) "geometrical" definition. Unfortunately, at the quantum level 
this operator displays a singular continuum limit, since it assigns non-zero topo- 
logical charge to lattice configurations extended over a size which scales with the 
lattice spacing (these configurations are called in the literature "dislocations" ) . 
An alternative approach g] consists in defining the lattice topological charge 
operator in terms of a local density, obtained from the naive discretization of 
the continuum operator. At the classical level, the resulting lattice operator, 
though not integer- valued, has the right continuum limit. At the quantum level, 
its matrix elements are connected to those of the continuum by non trivial renor- 
malizations |3| . This approach, relying only on first principles of Quantum Field 
Theory, is theoretically sound, but it has the inconvenience of the evaluation 
of the renormalizations, a theoretical uncertainty being introduced as a conse- 
quence in the lattice determinations. 

Among the cut off effects, those in the fermion sector are specially both- 
ersome. The Nielsen- Ninomiy a theorem states the impossibility of finding 
a local lattice action which is chiral invariant and describes only one species 
fermion in the continuum limit. A related problem is the impossibility of repro- 
ducing the right chiral anomaly in the continuum limit with chirally symmetric 
lattice actions. Since locality is of utmost importance for practical purposes, 
any lattice action for QCD must necessarily violate chiral symmetry. This vio- 
lation is so strong for the simplest discretizations of the fermion action ( "Wilson 
action" , "clover action" ) that all the chiral properties of the continuum Dirac 
operator are lost on the lattice, although they should be recovered in the contin- 
uum limit. An example of these properties is, at the classical level, the existence 
of zero modes with definite helicity, and their relation to the topological charge 
of the background gauge configuration through the Atiyah-Singer formula (In- 
dex Theorem) || : 
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Q = n L ~ n R 



(1) 



where Q is the topological charge and n L and n R are respectively the number 
of the left- and right-handed zero modes of the Dirac operator. 

At the quantum level, other important diseases are induced by the violation 
of chiral the symmetry; an example is the necessity of a fine tuning of the param- 
eters to recover the (spontaneously broken) chiral symmetry in the continuum 
limit (namely to get massless pions, chiral Ward identities, etc.). 

The freedom of choosing the lattice action among all possibilities in the 
universality class of the proper renormalized theory can be exploited in order to 
reduce the cut off effects at the price of rising the number of couplings involved 
in the action and their extension in space-time. The reward consists in more 
accurate physical results from simulations on coarse lattices. The drawback is 
of course the slowing down of the simulation procedure, due to the increased 
complexity of the action. Whether all this can be of practical or only academic 
relevance depends on the effective improvement of the accuracy of the results 
for a fixed cost in computational time. 

It follows as a consequence of the Wilson Renormalization Group (RG) ideas 
that there exist actions which give exact physical answers no matter how coarse 
the lattice is, the perfect actions They are located on the renormalized 
trajectory of any RG transformation. The first step toward perfect actions 
is represented by the so called fixed point (FP) actions, which are perfect at 
the classical level: this means, for asymptotically free theories, in the limit 
g — ► 0, where g is the (asymptotically free) coupling constant. It was shown 
by Hasenfratz and Niedermayer that the determination of the FP action of 
an asymptotically free theory is a feasible objective, since it requires only the 
solution of a classical minimization problem. The FP action is by construction 
perfect at the classical level, and as a consequence cut off effects are absent in 
the tree-level term of the perturbative expansion of the spectral quantities of the 
theory; for some specific observables, one-loop cut off effects are absent j|[ |9|, 
or strongly suppressed with respect to a standard discretization |jio| . 

In the context of the SU(3) Yang-Mills theory, the properties of the FP action 
are by now well known from a theoretical point of view and they have been 
extensively tested in Monte Carlo simulations [Jl2|, |l3|, |14|, with the result of 
very small cut off effects for spectral quantities even at moderate correlation 
lengths. Also, the ideas of the RG allow to build a well-defined topological charge 
operator on the lattice, the FP topological charge operator, which assumes 
integer values on all gauge configurations and is not affected by topological 
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defects @ |l6). 



The next step, the inclusion of quarks, is yet in a preliminary phase |l7|, p"8| . 
It is theoretically clear that many of the diseases of the simplest fermion 

actions produced by the violation of the chiral symmetry are cured by working 
with a FP action, whose (obligatory) chiral asymmetry is so mild that it has no 
effect in the physical results: the classical statements like the existence of zero 
modes of topological origin and the Index Theorem hold p3] ; at the quantum 
level, chiral symmetry is recovered at zero bare quark mass, without any tuning 
of an additional parameter, the currents get no renormalization, and there is no 
mixing between composite fermion operators corresponding to different chiral 
representations pifl . This properties are deduced using the Ginsparg- Wilson 
(GW) relation [|22|, which expresses the mildest breaking of the chiral symmetry 
for a "legal" fermion action; this relation is satisfied by any FP Dirac operator. 
Very recently, Liischer has pointed out that fermion actions verifying the GW 
relation have an exact continuous symmetry, which can be interpreted as a 
lattice implementation of the chiral symmetry without any contradiction with 
the Nielsen-Ninomiya theorem |2j| . 

It is worth to observe in this context that a solution of GW equation was 
obtained independently from the RG techniques. This solution is a product of 
the overlap formalism and shares with the FP Dirac operator its gentle chiral 
properties p4| . 

In this paper we want to verify the classical properties of the FP action in the 
fermion sector in a model simpler than QCD: the Schwinger model. The work 
is analytical for the pure gauge sector and numerical for the fermion sector. We 
investigate in particular how the approximation of the FP action with a finite- 
range fermion matrix affects its (classically perfect) properties. The paper is 
organized as follows. In Sec. 2 we recall the recursion relation obeyed by 
the FP Dirac operator, enlightening some properties of its classical solutions 
and bounds on the spectrum; in particular, we formulate the lattice version of 
the Index Theorem. In Sec. 3 we describe some topological properties of the 
continuum Schwinger model relevant for the subsequent work and define the 
RG transformation which generates the FP action under study, discussing the 
analytical results of the pure-gauge sector. Sec. 4 is devoted to the discussion 
of the numerical results concerning the spectrum of the FP Dirac operator. The 
paper ends with a summary of the conclusions in Sec. 5. 



2 The Fixed Point lattice Dirac operator 
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2.1 Recursion relations 



Let us recall schematically how the FP action for asymptotically free gauge 
theories can be computed; for further details, we refer the reader to the lit- 
erature [jl]] |6j. Consider a gauge theory, whose lattice action can be written 

as 



S[U,il>,#] = PS G (U) + $ n A nm (U)i> m , (2) 

where S G (U) is a lattice action for the pure gauge part, for example the Wilson 
single plaquette action, and A(£/) is a proper discretization of the continuum 
Dirac operator. On the lattice, the latter operator has the form of a matrix 
with space-time, Dirac and color indices, which in the above formula have been 
indicated in a collective notation. Since we do not want replica fermions in 
the spectrum, A(U) cannot be chiral symmetric, i.e. it must mix-up left- and 
right-handed spinor components. 

We define a physically equivalent action^] on a lattice with doubled lattice 
spacing through a RG transformation: 



exp [-f3'S G (U') - S'M'A'.U 1 )] = 

J [dUd^d^j] exp [-PS G (U) - ^A(77)V> - pK, G (U',U) 
-K F lCr(TP',TP',lP,lP,U)] , (3) 

where the primed fields are the degrees of freedom on the coarser lattice, defined 
through the gauge invariant kernels /C G and /C F . The parameters p and k f 
can be chosen arbitrarily^ in particular, we can take p = j3 K G> where k g is a 
fixed parameter. Note that in general the fermion action is not quadratic in the 
fermion fields after a RG transformation. 

In the class of theories under interest (including asymptotically free non- 
abelian gauge theories and the Schwinger model), the critical surface (where the 
continuum is attained) is at (3 — oo. The iteration of the RG transformation 
starting on this surface converges to a FP. When (3 — ► oo the integral on the 

2 With the same long distance physics. 

3 The choice is somehow restricted by the request that the RG transformation converges to 
a FP when iterated infinitely many times. 
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gauge degrees of freedom in the r.h.s. of Eq. (g) is saturated by the saddle point 
configuration; the solution of the recursion is then: 



S' G (U') = mm [S G (U) + K G JC G (U', U) ] , (4) 



S^',ip',U') = -In / [dtp dtp] exp [- VjA(L/ m '")V 



where f/ mm is the fine gauge field configuration which minimizes the r.h.s. of 
Eq. (|J). It depends on the coarse configuration U' . We have also (3' = (3 = oo. 
In this limit the problem is equivalent to a classical minimization problem plus 
a Grassmann integration. 

If the fermion kernel is quadratic in the fermion fields, the fermion action at 
f3 = oo remains quadratic after a RG transformation, as is evident from Eq. ([5]). 
We write in general 



where x b and x, y label the sites in the coarse and the fine lattices respectively; 
to is a matrix in space-time and color indices, depending on the gauge field U 
in such a way that gauge invariance is preserved, and normalized according to: 
^Zx^xtx = . The adjoint operation denoted with the symbol f acts on all 
indices, including space-time. 

The FP action is defined as 



5 FP = (3Sl P (U) + ^A FP (f/)V, (7) 

where S G r and A FP are the self- reproducing solutions of Eq. (Q) and Eq. (0) 
respectively. The recursion relation for the FP action of the pure gauge sector 
and for the FP Dirac operator read respectively: 



S pp ([/) = Sl P (U^(U)) + k g IC g (U,U—(U)) , (8) 



G 




(U) = Kf s l 



f u x b y b 




A FP + (t F wt w 



1 



') 



-t 



(9) 



where [/ min (?7) indicates the lattice gauge configuration defined on the fine lat- 
tice realizing the minimum in the r.h.s of Eq. (|J) for a fixed coarse lattice 
configuration U . Notice that, although not explicitly displayed, all the matrices 



The recursion relation for the inverse Dirac operator D = A -1 - the fermion 
propagator in the background of the gauge configuration U - is linear and there- 
fore simpler and often more convenientR: 



D ZA U ) = T,^( Um ' n ) D xl( U """)^yy b (U'-) + -S Xbyb . (10) 



2.2 Fixed point operators 

Suppose 0(U) to be some lattice discretization of a continuum operator 
depending only on the gauge fields. The RG transformation for such operator 



reads p5[ : 

XO'(U') exp [-/3'S' G (U') - S'M'^'iU 1 )] = 

[ [dU dip dip] 0{U) exp [- (3 S G {U) - ipA(U)ip - pK G (U',U) 



where the eigenvalue A is determined by dimensional considerations. This re- 
lation preserves (apart from a trivial multiplicative factor) the matrix elements 
of the operator 0(U) if, of course, also the action is changed accordingly. 

In the (3 — > oo limit, if the action is at the FP, the above equation reduces 

to: 

4 It should be noted however that, due to possible zero modes of the Dirac operator, its 
inverse does not always exist. 




K F JC F (l//,1p',lp,lp,U)] 



(11) 



7 



XO'(U') = 0(C/ m,n (C/')) 



(12) 



Eq. ( p"2| ) has a FP in the operator O fp (U) formally defined by the following 
limit: 

O fp {U) = lim (l/\) n O(U M (U)) , (13) 

n — >oo 

where U {n) {U) is obtained by applying n times - in a recursive way - Eq. (||) 
starting from the coarse configuration U: 

U m = U™ n (U ( - iy ) and U m = U. (14) 

As is evident from the definition (|l3|), the FP operator O fp is the solution 
of the FP equation: 

XO FP (U) = FP (C/ min (C/)) . (15) 

O fp (U) represents the classical perfect lattice operator, in the sense that its 
classical properties are the same of the corresponding continuum operator. 

A case of particular interest is that of the topological charge operator. In- 
deed, the FP charge operator Q FP {U) allows a theoretically sound definition of 
the topology on the lattice, since it assumes integer values and it is not affected 
by short-range fluctuations ( "dislocations" ) |l5|, [l(| . 

2.3 Classical solutions and the Index Theorem 

In analogy with the pure-gauge sector, where there exist exact classical so- 
lutions of the lattice equations of motion having the main properties of the 
continuum solutions (e.g. scale invariance) it is also possible to prove jl9| 
the existence of exact solutions of the lattice FP Dirac equation, given by the 
zero modes of the FP Dirac operator. The arguments given here follow p9[ ; see 
also [^0| for an alternative approach. 

Suppose f<£} to be a solution of the FP Dirac equation in the background of 
some coarse configuration U (c) (Dirac and color indices are understood): 
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£A^(C^)/<°> = . (16) 

Vb 

Using ([)[) it is easy to verify that 



^ 1 ] 4« B (u-)m (it) 



yy B \ \~ \~ / \~ I / Xy 



A FP (C/™) + K F wt(U-) w ({/™) 



is also a solution of the Dirac equation in the fine lattice with the background 
configuration C/ min ([/ (c) ): 



E 



AH(U-(U^)) f y f > = (18) 



Conversely, let f!j. f) be a solution of the Dirac equation in the fine lattice, 
with the background configuration C/ mi "([/ <c) ). Then 



m = £^([/™)^> (19) 



is a solution of the Dirac equation in the coarse lattice with the background 
configuration U {c) . 

Eqs. dl7| ) and (|l^) establish a one-to-one correspondence between the zero 
modes of the FP Dirac operator in the background of the coarse gauge configura- 
tion U [c) with those of the same operator, but considered now in the background 
of the fine gauge configuration C/ min ([/ (c) ). From Eq. ( [l9| ) we see that two re- 
lated zero modes must have the same chirality, since ui does not depend on 
Dirac indices. By iteration, one can relate in the same way the zero modes of 
A^([/ (c) ) to those of AH(m n) (U^)), where f7 Cn) (E/ (c) ) is defined as in Eq. (|TJ). 
In the limit of infinite iterations n — > oo, f7 (n '(t/ (o) ) converges to a continuous 
configuration. The topological charge of such configuration is given by 



Q = lim L (C/ ( " ) ([/ (c) )) , (20) 
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where Q L indicates any proper lattice discretization of the topological charge 
operator. The r.h.s. of the above equation is just the definition of the FP topo- 
logical charge of the coarse lattice configuration U (c) , Q FF (U (c) ) (see Eq. (|l3|); 
in the case of the topological charge, a dimensionless operator, A = 1). For 
the limiting continuous configuration the Index Theorem holds, since A^([/ <n) ) 
converges to the continuum operator when n — > oo. Using the above estab- 
lished one-to-one correspondence between the zero modes, and the fact that 
this correspondence preserves chirality, we arrive to the conclusion that (l9[ |20| : 

i) the FP Dirac operator A FP (£7) possesses exact zero modes of topological 
origin; 

ii) these zero modes have definite chirality and satisfy the Index Theorem, if 
the lattice topological charge is the FP topological charge Q Fr (U) of the 
background lattice gauge configuration U. 



2.4 The Ginsparg-Wilson relation 

Let us show now how, forced by the Nielsen-Ninomiya theorem, the FP Dirac 
operator violates the chiral symmetry (n], |2(|. Consider the n-th iteration of 
Eq. ©: 

D HvS u ) = E n &(tf) Dll(u ln) ) n&Hu) + <y b (U) , (21) 

xy 



where 



= £ + E E «s- n^(co) (22) 



and 



«Sx(tO = E ^Aun^ lX2 (u (2) )...^ 1 Au ( * ) ), (23) 

Xl . ..Xi-l 

with E/«(f7) defined as in Eq. (111). 
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The background configuration U {n) in Eq. ( plf ) becomes arbitrarily smooth 
when increasing n. Therefore, when n — > oo the propagator in the r.h.s. of ( pl| ) 
converges to the continuum propagator, which is chiral invariant. Since f2 ( ™' 
commutes with 75 (we recall that u> does not act on Dirac indices), we have |L9| 



{^,75} = 2 75J R- b (C/), (24) 
where i? FP the limit n — > 00 of R ln) . Equivalently, we have 

{A- , 75 } = 2 ( A FP 75 R ¥F A pp ) . (25) 



The hermitian matrix i?^ p yb is local and its spectrum is bounded 



1 < K J^1< J_ , <«<l; (26) 
(u,u) 1-g 

the lower bound comes trivially from Eq. (|22]), while the upper bound (specif- 
ically the constant q) is RG-transformation-dependent. In the case of the 
overlapping-symmetric fermion RG transformation |l7| |2l| , and its gauge-invariant 
extension for the interacting theory B here considered (see the following), one 
has q = 1/2. 

Ginsparg and Wilson proved [p2| that any lattice discretization of the Dirac 
operator which verifies Eqs. of the form (p4])-(p5|) - it is crucial the locality of 
the matrix RZ P , , - reproduces the correct axial anomaly in the continuum limit. 
We refer Eqs. (|24D-(|25|) as the Ginsparg- Wilson relation (GWR). It expresses 
the explicit violation of chiral symmetry by the lattice regulator. The fact that 
this violation takes the form of the GWR ensures that many important chiral 
properties of the continuum Dirac operator hold also for its FP discretization. 

For instance, using the GWR the Index Theorem has been proven in Ref. 
directly on the lattice, without any reference to the Theorem of the continuum, 
finding also the following relation between the FP topological charge and the 
FP Dirac operator: 



\U) = -itr( 75 {A pp ([/), J R pp (L/)}). (27) 
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2.5 About the spectrum 

The GWR strongly constrains the spectrum of the FP Dirac operator. Con- 
sider first the zero modes. Let us denote by v a vector column which contains 
the spatial, Dirac and color indices. 

Property: 

v is a zero mode of A FP if and only if 75 v is. Therefore, the zero modes can be 
chosen eigenstates of 75. 

The proof is trivial using ( p5| ) : 

A FP l5 v = -7 5 A FP i; + { A FP , l5 }v = . (28) 

In what follows we shall prove that the FP Dirac operator is bounded, and 
we shall find the explicit bound on the spectrum. First, we notice that the FP 
Dirac operator verifies the hermiticity condition 

A Fp t = 7 5 A FP 75 , (29) 

since it is preserved by the recursion relation (^). The hermiticity property, 
which also holds for the Wilson Dirac operator, implies that the eigenvalues of 
A FP are either real or comes by pairs of complex conjugate numbers. It also 
implies that if v is an eigenvector corresponding to a non real eigenvalue A 
(A =^ A*) then 

faw) = . (30) 
Multiplying on the left by 75 both sides of ( p5| ) and using (]29|), we obtain: 

A fp + A FP t = 2A FPt i? FP A FP . (31) 

Let v be an eigenvector of A FP with eigenvalue A, normalized to unity: (v, v) = 1. 
Multiplying both sides of ([H]) by t>t on the left and by v on the right, we find: 

A + A* = 2(v,R FP v) |A| 2 . (32) 
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The above equation, in conjunction with the bounds ( pq ) for the R FF matrix, 
allows to state analogous bounds for the spectrum of the FP Dirac operator p(J] : 

-|A| 2 < A + A* < 2 |A| 2 . (33) 
K F (1 — q) K F 

These inequalities mean that the spectrum of the FP Dirac operator lies in the 
complex plane inside a circle of radius k f /2 centered at (k f /2, 0), and outside 
a circle of radius (1 — q) K F /2 centered at ((1 — q) K F /2, 0). 



3 The case of the Schwinger model 



3.1 Topological properties of Schwinger model 



The Schwinger model is a good laboratory to test all these ideas, since 
it is much simpler than four-dimensional non-abelian gauge theories, but still 
its gauge sector has non-trivial topology. The topological charge of a gauge 
configuration is given by: 



Q = ± I d 2 xF 12 , (34) 

where: F 12 = d 1 A 2 — d 2 A 1 and e is the electric charge. For the theory defined on 
a torus T 2 ( x^ S [0, L J, /i = 1,2), topologically non-trivial configurations must 
necessarily have a jump at the boundary. Of course, gauge invariant functionals 
of the gauge fields must be continuous. The simplest topologically non-trivial 
gauge configurations are solutions of the classical equations of motion, therefore 
with constant F 12 ; in this case a possible expression for the gauge fields is |27| : 



Af at (x) = -$x 2 , Af at (x) =0, F 12 = (35) 

The fields A™ st (x 2 = L 2 ) and A i { 1 '* t (x 2 = 0) can be connected by a transition 
function (ensuring continuity for gauge-invariant quantities) if $ — 27rQ/(LiL2e), 
with Q some integer number, corresponding to the topological charge of the 
gauge configuration on the torus according to definition (|34|). 

In perfect analogy to the case of QCD, also in the Schwinger model it is 
possible to relate the chiral properties of the zero-modes of the massless Dirac 
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operator to the topological properties of the background gauge configuration 
A^. This relation is given by the Atiyah-Singer Index theorem, Eq. (Q). A 
peculiarity of two space-time dimensions is that: 



Q{A) > => n R = , 

Q(A) < n L = , (36) 

the so called "Vanishing Theorem" (^8|. So, in the case Q(A) ^ 0, the zero 
modes have all the same chirality, given by the sign of Q(A). 



3.2 The pure gauge sector on the lattice 

To be able to study the topological properties of the Schwinger model on the 
lattice, we discretize the theory in the compact formulation, where the lattice 
gauge fields are described by angular variables a^{x) G (— it, tt] and, forced by 
the necessity of preserving gauge invariance, the fermions are coupled to the 
gauge fields through U{\) elements of the form: 

U^x) = e ia »W . (37) 

The lattice action must be a periodic function of a M . Within this formulation, 
we can implement periodic boundary conditions also for the topologically non- 
trivial gauge field configurations. The lattice discretization of the continuum 
instanton solution on the torus ( |35| ) is p9[ : 



af et (x) = -$ L x 2 (mod27r) 

a™\x) = if x 2 = 0,l,...,AT 2 -2 

= $'JV 2 x 1 (mod 2tt) if x 2 = N 2 - 1 , (38) 

where $ L = 2ttQ/{N 1 N 2 ) and [9 (mod27r)] G (-vr, n}. 

To carry on our program of finding a FP lattice action we must define a RG 
transformation. For the pure gauge part, we consider the following kernel (the 
primed variables are those relative to the coarse lattice): 
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K, G (a',a) = ^2 ([ a 'iJ,( x b) - ^n( x b)] (mod27r)) 2 , (39) 



where r^(x&) — [a (1 (i) + a^x + ft) ](mod 2tt) . 

In the limit k g — > oo, Eq. (^) reduces to a problem of constrained minimum: 

S' G (a') = 4 min {S G (a) \ T^x b ) = a'^x b )} . (40) 

{a(x)} 

The 4 factor comes out since in two dimensions the electric charge is a relevant 
parameter having the dimension of a mass and, therefore, the coupling constant 
is trivially renormalized at the lowest (tree level) order: /?' = /3/4. 

In the Appendix we proof that the FP of the above recursion relation is 
given by: 

So» - 5 E[ F -(z) (mod2^)f , (41) 

X 

where F^ v (x) is the lattice field strength tensor: 

Fjlu = a n( x ) + oc v {x + p) - a^(x + v) - a u (x) . (42) 

We have so obtained the Manton action |30| of the U(l) theory. This is 
not a surprise since it is known pH that the U(l) pure gauge theory in two 
dimensions is equivalent to the one dimensional quantum rotor, for which the 
Manton action was proven |3^] to be classically perfect. In |53| it was pointed 
out that the string tension computed with the Manton action shows perfect 
scaling. 

In the Appendix we also report an explicit solution (among all the possi- 
ble gauge-equivalent ones) for the fine configuration minimizing the r.h.s. of 
Eq. ©, a min (c/), giving trivially through @ C/ mi " ([/'). 

The action in Eq. ( ^l| ) can be rewritten in terms of the link variables: 
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S l"( U ) = -£(Imln(C/ 12 (aO)) 2 

X 



(43) 



where U 12 (x) represents the usual product of links around the 1 — 2 plaquette. 
The r.h.s. of (^) is not defined for the so-called exceptional configurations^], 
which have U 12 (x) = — 1 for some x. These, however, have zero measure in the 
path integral and can be ignored. We see that the FP action, though ultra-local, 
involving only first neighbor interactions, has a non-polynomial dependence on 
the link variables. In particular, considering the expansion: 



ln(C/ 12 ) = -£r(l-£/ 12 ) ft , (44) 
fe=i 

we understand that the fixed point action takes contribution from monomials 
of the link variables of arbitrary order, with a slow suppression with the in- 
creasing order. This may suggest that, in the framework of the non-abelian 
Yang-Mills theory, the choice to parametrize the FP action by a finite set of 
operators of loops of the link variables is definitely not the most suitable, and 
other approaches should be dared jl9| . 

The lattice configuration of Eq. (^8|) is a classical solution for the FP ac- 
tion (jFL|)-(|43|). It is in general expected that a FP action (which - we recall 
- is a classically perfect action) reproduces the continuum action on classical 
solutions Q. This is indeed the case in the present context, since, as it is clear 
from the expression (|4l|): 

PS? (ff*) = 2 -^- ; (45) 

exploiting the relation = l/e 2 a 2 , the r.h.s of the above equation can be 
rewritten in physical units as 2ir 2 Q 2 / L 2 e 2 , i.e. the continuum action of the 
configuration (|35|). 

To end the discussion of the pure gauge sector, we shall study the FP topo- 
logical charge. Using the expression given in the Appendix for a min (a'), it is 
straightforward to prove that the solution of the recursion relation: 

5 The term exceptional configurations is used in a different context to denote those con- 
figurations for which the inversion of the Dirac operator is numerically complicated by the 
presence of nearly zero modes. 
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Q F » = Q FP (a-""(a)) (46) 

is given by (cfr. p2|): 



2 P » = ^^^W(mod27r) . (47) 



The r.h.s. of the above equation is an integer, since Ff 2 (x) = due to the 
periodic boundary conditions. Notice that the geometrical charge, 



Q° = ^ E Im M^) - ( 48 ) 



coincides with the FP charge for non-exceptional configurations^ Since excep- 
tional configurations have zero measure in the path integral, we see that the 
geometrical charge gives the correct continuum limit. This is somehow different 
from what happens in the non-abelian theories, where the geometrical charge is 
singular in the continuum limit due to the dislocations p3|. 



3.3 RG TRANSFORMATION FOR THE FERMIONS 

For the fermion sector, we choose the quadratic RG kernel (0) with the 
following matrix ui H : 



\[2 { 1 

UJ Xb x(U) = —<S x ,2x b + 2 E U »( 2x b) 4,2x 6 +A + 
I ju=±l,±2 

[U t ,{2x b )U l/ (2x b + fi) + U u (2x b )U t ,(2x b + i))}5 x ^x b +f i +A 
where, as usual, U_i(x) = U'(x — 1) and the analogous relation for 1 — > 2. 

6 The FP charge fixes the prescription: ln( — 1) = in in the case of an exceptional 
configuration. 



(49) 
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Once the gauge sector has been worked out, with a prescription for £/ min (£/), 
we are in a position to solve, for a given background lattice configuration U, 
the iteration equation for the FP Dirac operator, Eq. @. Of course, due to 
the complexity of the equation, this cannot be done analytically, and so we 
have to turn to some approximate method. In the next Section we describe 
the technical details of the numerical method and discuss the reliability of the 
approximations . 



4 The numerical computation 



The ideas of Sec. || would be of only academic relevance if we were not able to 
compute an approximation of the FP Dirac operator A FP in terms of a matrix 
with a finite (and sufficiently short) interaction range. This approximation 
should be reliable, in the sense that all the theoretical properties mentioned in 
Sec. U should be preserved with good accuracy. We shall show in the remaining 
of the paper how this can be done in practice. 

We have computed the FP Dirac operator (or better, some ultra-local ap- 
proximation to it, see below) for several background gauge configurations. We 



have considered the instanton configurations of Eq. (38) for Q = 1 and Q = 2: 
we will refer to them as to smooth configurations. We have also studied the 
effect of the super-imposition on these smooth configurations of a random per- 
turbations with a given size. This has been done by adding to the gauge field of 
the smooth configuration a perturbation of the form 5A^ = size x r, where r is a 
random number between —1 and 1. Finally, we have also considered completely 
random configurations. 

We solve the FP equation (0) by iterations. Our strategy is the following: 
we start from some lattice definition of the Dirac operator and, for a fixed 
background configuration U, we let it evolve under RG according to Eq. (^|). 
After a large enough number of RG transformations, the running Dirac operator 
will converge with sufficient accuracy to the FP Dirac operator. 

In order to obtain, according to Eq. (^|) , the Dirac operator evolved through 
n RG transformations A {n) (U), defined on a N x N toroidal lattice, we proceed 
in the following way: first of all we construct by iterative procedure the se- 
quence of gauge field configurations: U {0) {U) = U, U {1) (U), . . . , U (n) (U), living 
respectively on a N x N, 2N x 2N, . . . , 2 n N x 2 n N toroidal lattice. Eq. @ 
allows us to calculate A (1) (U^'^QT)), defined on a 2 n ~ 1 N x 2 n-1 JV lattice, by 
putting A (0) ([/<"' ([/)), defined on a 2 n Nx2 n N lattice, on the r.h.s. of the equa- 
tion; once obtained A (1) (C/ ( " _1) (£/)), the procedure can be iterated to determine 
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A (2) (C/ ( " 2) (U)) ... (defined on smaller and smaller lattices) up to A {n) (U) on 
the final N x N lattice. 

Considering that the dimension of the matrices entering the above described 
procedure explodes when increasing n (taking into account both space-time 
and Dirac indices, the largest matrix involved, A m (U (rl} ({/)), has dimension 
2 2n+1 N 2 ) it is clear that some approximation must be introducecQ. In what 
follows we discuss these approximations and the convergence of the iterative 
procedure. 



4.1 Approximations 

From RG theoretical arguments, we expect the FP Dirac operator to be 
local, though coupling fermions at arbitrarily large distances. In this context 
locality means that the matrix elements are exponentially suppressed with 
the distance \x — y\. We know from a perturbative computation |J that the 
most important couplings are those inside the region \x /J , — y l j,\ < 2. This local 
structure is very important from the practical point of view, since it allows us 
to escape the memory-space problems connected with the allocation of matrix 
elements of the Dirac operator on the finer lattices. 

Therefore, as an approximation we consider only ultra-local fermion matri- 
ces: we neglect in all stages of the recursive procedure the matrix elements 
which couple fermions at a distance larger than a certain maximum range r g . 
In other words, we make a sharp cut in the interaction, considering a fermion 
action of the form: 

x b x b +l i>x„+i , (50) 

x b , I 

where l^, fi = 1,2, takes values between —r g and r g . It is worthwile to stress 
that this ultra-local approximation maintains all the symmetries of the exact 
action, including hermiticity (Eq.( |2^)). To study the locality of the FP Dirac 
operator and to keep the effects of the cut under control, we consider the two 
cases r g = 1 and r g — 2. 

Even though within the ultra-local approximation memory problems are 
partially solved, we have still to cope with the CPU time, which grows with the 
number of iterations as ~ 2". As a consequence, we are not able to make many 

7 Note that in Eq. M a matrix inversion is required. 
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iterations, and in practice we have to stop at n = 4 or n = 5. This means that 
we do not get the exact FP Dirac operator, but some asymptotic approximation 
to it. 

At the level of free fermions, the sharp cut of the couplings has the effect, 
besides an unphysical change of the normalization of the fermion fields, of gen- 
crating a fermion mass, with the consequence of bringing the action outside the 
critical surface. Being the mass a relevant parameter, the convergence towards 
the critical FP is therefore spoiled. 

Explicitly, writing the free Dirac operator as 

A ?y = J2 p^ x - y) > + A ( s - y) > ( 51 ) 

we have the following normalization conditions: 

X X 

the first equation sets the normalization of the fermion field, and the second 
is the massless condition. To preserve conditions (|5^) in the course of our 
approximated iterative procedure, we must add to each step a multiplicative 
renormalization of the fields and an additive renormalization of the mass. 

In the interacting case, we use the renormalization constants computed in 
the free case. This ensures that the iterative procedure keeps the fermion action 
in the universality class of the massless Schwinger model. 

The renormalization constants are a measure of the goodness of the ultra- 
local approximation: the closer are the field and mass renormalization respec- 
tively to one and zero, the better the approximation. Table displays the 
rcnormalizations as a function of the number of iterations for the cases under 
study r g = 1 and r g — 2. We see that for r g = 2 the renormalizations are 
very small ( O(10 -3 ) ), as expected from the locality properties of the FP Dirac 
operator. 
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r 9 


= 1 


T 9 


= 2 


step 


field 


mass 


field 


mass 


1 


1.1957 


-0.7431 


1.0026 


0.0038 


2 


1.1661 


-0.7467 


0.9998 


0.0040 


3 


1.1508 


-0.7477 


0.9999 


0.0022 


4 


1.1428 


-0.7482 


1.0019 


0.0002 


■5 


1.1381 


-0.7491 


1.0018 


-0.0020 



Table 1: The values of the field and mass renormalization constants after each 
RG step, for the two sharp cuts considered. 



4.2 Convergence 

We call A (0) (C7) the starting discretization of the Dirac operator. In order 
to improve the convergence of the RG iteration, it is important to choose such 
operator as close as possible to the FP. For the construction of an approximate 
expression of the FP Dirac operator, the analytical knowledge of the FP matrix 
for the free fermions, A^ froB , has been exploited. The idea is to transform the 
couplings between two free fields, ip x and ip v , into gauge-covariant couplings 
between fermions interacting with a gauge field; the simplest way of doing this 
consists in multiplying A^ tro ° by the lattice parallel transporter on the shortest 
path^l from x to y. The result is our definition of A m (U) xy . In order to avoid 
complicated paths, we have restricted the couplings of A m (U) x y inside the 
region |a? p - scJJ < 1. 

As a criterion to study the convergence of the iterative process, we have 
chosen the operator norm: 

IWU =-p{¥^} . (53) 
M I IMI2 J 

where v is any vector and ||w||2 = \/Yl v i v i- ^ n Table |2| we show the values 
of the norm of the difference of two Dirac operators, corresponding to two 
consecutive iterations, in the free case and, in the interacting case, for the 
smooth background configuration given by Eq. (pq) , random fluctuations around 
it (sizes 0.5 and 1), and a completely random configuration; the background 
configurations have all fixed point topological charge Q FP = 1. The lattice size 
is N = 6 (except for the case of the random configuration, where N = 4) and the 

8 In case of more than one path, an average that keeps the lattice symmetries has been 
taken. 
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i 


free 


smooth 


pert. 0.5 


pert. 1 


random 


1 


0.3320 


0.3248 


0.4138 


0.5989 


1.9678 


2 


0.1997 


0.1971 


0.2068 


0.2276 


1.6476 


3 


0.1148 


0.1143 


0.1184 


0.1284 


1.4608 


4 


0.0630 


0.0631 


0.0661 


0.0740 


0.7542 


5 


0.0348 


0.0352 


0.0382 


0.0461 


0.4112 



Table 2: The value of || A (i) — A (i-1) ||oo for two consecutive iterations, i — 1 and 
i (i = is the starting operator). The case considered is N = 6, r g = 1 and 
qff _ ^ p Qr ranc | om configuration, N = 4. 

range of the interaction r g = 1. In all the other cases (different lattice size, range 
r g and topological charge) the results are qualitatively similar. Quantitatively, 
the deviations from the numbers of Table ^ are not large, ~ 10% at most. 
Observe that the data of Table || indicate, as expected, the presence of irrelevant 
operators of dimension three. 

We have observed that the rougher the configuration is, the worst the con- 
vergence (the numbers of Table ^ are an example) . This is natural, since for a 
rough configuration U more iterations are required in order that the configura- 
tion at the finest level, U in) (U), becomes smooth enough to ensure convergence. 
This explains in particular the quite slow convergence for completely random 
configurations, as is manifest in Table 0. 



5 Numerical results 

In this Section we discuss the numerical results coming from the checks of 
the classical properties of the FP Dirac operator. 

In the iteration of the FP equation a matrix inversion must be performed 
at each step. We use for this an recursive algorithm, which gives results with a 
precision of 10 -4 . This means that in the following, numbers smaller 10~ 4 are 
compatible with zero. 



5.1 Spectrum and bounds 

We start by discussing the properties of the spectrum. We use k f — 4, since 
in the free case this value gives the most local action. Eq. (|3^) for q = 1/2 tells 
us that the spectrum must be inside the circle of the complex plane of radius 
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2 centered at (2,0) and outside the circle of radius 1 centered at (1,0). While 
in no case we have found a violation of the first bound (it is true for each value 
of the number of iterations, in fact, it turns to be true also for the starting 
operator A (0) ), the second one holds with good approximation only for r g = 2. 

Fig. [j] displays the spectrum of the FP Dirac operator on a 6 2 lattice for 
background configurations with Q FP = 1. The Dirac operator has been obtained 
through 5 RG iterations. The deformation of the spectrum induced by the 
cut can be observed by comparing Figs. [j]-a and Q-b, which refer to the two 
different cuts r g — 2 and r g = 1 respectively, in the case of a smooth background 
configuration. We see that for r g = 2 the bounds on the eigenvalues are well 
satisfied, while for r g = 1 a slight violation comes into play in the left edge 
of the spectrum due to the approximation of the sharp cut of the couplings. 
Comparison between Figs. 0-b and 0-c,d allows to evaluate the stability of the 
spectrum against perturbation of the background configuration, corresponding 
the two latter cases to a perturbation of the background configuration of size 
0.5 and 1 respectively (the cut is always r g = 1). In all cases, except d, we find 
three real eigenvalues A < Ai < A 2 : the lowest one, A corresponds to the would 
be zero mode, see below; X 1 is twofold degenerate. In the case d, where a strong 
perturbation of the background configuration is applied, the real eigenvalues are 
six, with no degeneracy. The higher real eigenvalues, lying on the right edge 
of the spectrum, are always decoupled from the infrared part of the spectrum, 
where the lowest real eigenvalue A is located. 

It is interesting the comparison with the case of the (massless) Wilson oper- 
ator, which has been studied recently in Ref. |35|; the eigenvalues fulfil in this 
case an additional symmetry p5[ , the reflection with respect to the point (2,0) 
in the complex plane^. For comparison with the FP case, we report in Fig. g the 
spectrum of the Wilson operator for the smooth configuration with topological 
charge one. 

In Fig. H the evolution of the spectrum under the RG process starting from 
the Wilson operator is displayed, for n — 1, 2, 3 and 4, in the case of the smooth 
configuration. We see that already for n = 1 the spectrum has assumed the 
main features of its characteristic geometry, and only small changes take place 
thereafter; in particular, for n = 4 the convergence to the final spectrum seems 
to be reached with good accuracy. It is also possible to observe the progressive 
decoupling of the two higher real eigenvalues Ai and A 2 , drifting towards higher 
values, from the lowest one A , converging to zero. 



9 Actually, this additional symmetry holds only for even lattices. 
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5.2 Zero modes and index theorem 



From the discussion of Sec. g we know that the FP Dirac operator must 
have zero modes in the topologically non trivial sectors, with definite chirality 
and satisfying the Index Theorem, Eq. (Q) , if the continuum topological charge 
is replaced by the FP topological charge Q FP . In the Schwinger model the 
Vanishing Theorem ( |36| ) must also be true. As a consequence, in the case 
Q FP > 0, the FP Dirac operator should have Q FP zero modes with positive 
chirality. Since we do not attain exactly the FP, we do not expect exact zero 
modes. Rather, we expect that the Q FP lowest lying eigenvalues converge to zero 
when the number of iterations tends to infinity. Still, the cut of the couplings 
in the region \x' — x\ > r g can produce unwanted effects. 

We start the discussion with the smooth configuration. In Table || we can see 
the eigenvalues of smallest modulo as a function of the number of iterations i. 
We report also the overall chirality of the zero modes (trace of 75 in the subspace 
of zero modes, X^( w '75 w ))' which should equal the FP topological charge of 
the background configuration. Notice that the eigenvalues are (within machine 
precision) real for i = 0; this is in agreement with the observation J35| that 
eigenvectors with non-zero average chirality must have necessarily real eigen- 
values. This is an exact statement, relying only on the hermiticity property of 
the fermion matrix (|2^), which is preserved by the RG procedure. The small 
numbers O(10~ 6 ) found for the imaginary part of the eigenvalues for i > 1 
are due to the approximation introduced when operating the matrix inversion 
in the recursion (the ultra-local approximation does not violate the hermitic- 
ity property). In accordance with the Index and Vanishing Theorem, in the 
case Q FP = 1 we find only one non-degenerate small real eigenvalue, while for 
Q FP = 2 the lowest lying eigenvalue is twofold degenerate. The real part of the 
eigenvalues converges towards zero when increasing the number of iterations; 
the law is roughly as 1/2*, as expected from the dimensional analysis of the 
lowest dimensional irrelevant operators. This law breaks down for i = 5 in the 
case r g = 1, where we find numbers smaller than expected. The eigenvectors 
have in any case definite positive chirality, regardless to the number of itera- 
tions. In fact, this is already true for the starting Dirac operator, and the cut 
of the couplings seems not to be effective in this regard. 

In conclusion, our results are consistent with the picture of a FP Dirac 
operator exactly verifying the Index Theorem, and being well approximated by 
ultra-local fermion matrices. 

It is interesting to compare these results with those given by the Wilson oper- 
ator, displayed in Table ||. Notice how the Index theorem is only asymptotically 
verified for N — > 00, i.e., in the continuum limit. 
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The numbers of Table |5| have been obtained by adding a random pertur- 
bation of size 0.5 to the smooth configurations previously considered (see also 
Table |6|) . Notice that the lowest lying eigenvalues are about twice larger than 
the corresponding for the smooth configuration. Again they decrease by a fac- 
tor of about 2 after each iteration. Notice also that the random perturbation 
removes the degeneracy of the two lowest lying eigenvalues in the case Q FP = 2. 
Finally, we see that the chirality is no more perfect from the start, but already 
for i — 4 it is compatible with 1 within our numerical accuracy. 



5.3 The topological charge 



The FP topological charge can be obtained directly from the FP Dirac oper- 
ator through the formula (|2^). This relation could be useful for the calculations 
of topological quantities in unquenched QCD with the FP action, since it gives 
automatically a prescription for the FP topological charge, once a parametriza- 
tion of the FP fermion matrix has been found. In our approach, we can check 
how accurately this relation is verified as the number of iterations of the FP 
equation is increased. The results are displayed in Table 0. We see that the 
difference with the FP topological charge is less than 14% after 4 iterations and 
less than 4% after 5 iterations with r g = 1 in the case Q = 1. Moreover, the 
results are stable under the perturbation of the background configuration. With 
r g = 2 we observed only a slight improvement. 



6 Conclusions 



We have presented an explicit example which illustrates that it is possible 
to compute an ultra-local approximation of the FP Dirac operator which pre- 
serves up to a very good accuracy all the important properties of the exact 
FP operator. We showed by numerical computation that in the case of the 
Schwinger model it is a very good approximation to consider only couplings 
between fermions contained in the same size-two plaquette (translated in 4 di- 
mensions: size- two hypercube). We have confined our work to the classical 
properties, like zero modes and their chirality, topological charge and the Index 
Theorem. We have checked these properties in the case of instanton solutions 
of the equation of motion and random perturbations around them. We believe, 
however, that the quantum properties of the FP action, like for example the 
non-renormalization of the pion mass, will also hold with similar accuracy. An 
indication of this is that another property derived from the GWR, namely the 
relation (27) between the FP Dirac operator and the topological charge of the 
background gauge configuration, holds within this same accuracy also when in- 
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troducing quantum fluctuations around the classical instanton solutions. We 
are nevertheless conscious of the fact Q that these mildly perturbed configu- 
rations are not representative of the whole MC set of configurations at thermal 
equilibrium. 

Recently a great attention has been paid to the Index theorem on the lat- 
tice. The goal was to find some traces of it, cither using the standard Wilson 
action |55|, ||(| , some chiral improved version of it |57| or approximated solutions 
of the GW relation constructed from the overlap formalism f3S[| . In all these 
cases, the Index theorem seems to hold only in a probabilistic sense. We hope to 
have convinced the reader that the exact Index theorem of Ref. holds very 
accurately even with an ultra-local approximation of the FP Dirac operator. 

Our ultimate conclusion is therefore that a parametrization of the FP action, 
involving only very localized couplings, with good scaling properties should be 
possible. Moreover, following the suggestion by Hasenfratz and Niedermayer, 
one can use the freedom typically present in this kind of fits to enforce up to 
a very high precision the classical properties for some instanton-like smooth 
configurations. A parametrization of the FP action of the Schwinger model was 
already found in |39| , but locality was not the first priority of the authors in that 
work. Really, this is urgent problem only in view of four dimensional theories. 

We solved the FP equation by iterations. The analytical solution of the 
pure gauge part (this is an extra-bonus of the two dimensions) enabled us to 
concentrate the numerical effort in the fermion sector. Also the low dimen- 
sionality of the model was of great help, but still, it was not possible to make 
many iterations. With five we got a good precision. In four dimensions, how- 
ever, memory-space and CPU-time problems are likely to be an obstacle to this 
naive iterative approach. Different strategies to solve the FP equation should 
be devised. A possibility to avoid extremely large lattices is to use, as in [^3, a 
parametrization of the FP action from the very beginning, namely to solve the 
iteration equation in a finite space of couplings. 
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Table 3: The lowest lying eigenvalue and the chirality of the corresponding 
eigenvectors of the Dirac operator for a smooth background gauge configuration 
(Q = 1 and 2), for each number of iterations i. The lattice size is iV = 6. r g is 
the range over which the approximated action extends, see subsection 1.1. 
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Table 4: The lowest lying eigenvalue and the chirality of the corresponding 
eigenvector of the Wilson Dirac operator on lattices of different sizes N, with 
the smooth background configuration having Q = 1. 
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Figure 1: The spectrum of the FP Dirac operator A FP . The lattice size is 
N = 6 and the number of iterations n = 5; a) smooth background configuration 
with Q FF = 1; the range of the couplings considered in the action is r g = 2 



(see subsection 4.1); the two circles represent the bounds on the spectrum, 



as explained in the text; b) as in a), but with r g = 1; c) and d) r g =l and 
background configuration with perturbation of size 0.5 and 1 respectively (still 
Q FP = l). 
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Figure 2: The spectrum of the Wilson operator A w on a 6 2 lattice for the 
smooth background configuration with Q FP = 1 . 
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Figure 3: The evolution under RG of the spectrum of the Dirac operator on a 
6 2 lattice starting from the Wilson discretization A w , for number of iterations 
n = 1 (a), n = 2 (b), n = 3 (c), n = 4 (d). The smooth background configuration 
has charge Q FP = 1 and r g = 1. 
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10" 


-6 


0.9998 


4 


0.0299 


7.9 


io- 


-6 


0.9998 


0.0174 


-1.5 


io- 


-7 


1.0000 




0.0490 


-4.4 


10" 


-6 


0.9992 


0.0133 


3.5 


10" 


-6 


1.0000 



Table 5: The same as table [3| adding a random fluctuation of size 0.5 to the 
smooth configuration. 



i 


Re A 


ImA 




(V, 15V) 





0.3348 


2.0 10- 


■16 


0.9282 


1 


0.1824 


1.1 io- 


-5 


0.9869 


2 


0.1245 


9.7 10" 


-6 


0.9963 


3 


0.0908 


5.7 10- 


-6 


0.9986 


4 


0.0679 


8.7 10- 


-6 


0.9990 


■5 


0.0520 


-2.2 10' 


-6 


0.9989 



Table 6: The same as table y adding a random fluctuation of size 1 to the 
smooth configuration. In this case Q = 1 and r g = 1. 
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= 1 



= 2 



i 


smooth 


Pert. 0.5 


pert. 1 


smooth 


Pert. 0.5 


Pert. 1 


1 


0.2289 


0.2280 


0.2118 


0.4565 


0.4498 


0.4263 


2 


0.5114 


0.5067 


0.4906 


1.0067 


0.9954 


0.9691 


3 


0.7223 


0.7185 


0.7100 


1.4143 


1.4020 


1.3902 


4 


0.8668 


0.8637 


0.8663 


1.6959 


1.6839 


1.6878 


5 


0.9631 


0.9608 


0.9740 


1.8844 







Table 7: The value of — t^ - tr (7 5 {A FP , i? FP }) for the smooth eonfigurations 
and perturbations around them; N = 6 and r g = 1. 
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A Appendix 



In this Appendix we solve explicitly the minimization problem ( p0| ) when 
in the r.h.s. the Manton action is taken. The result will be a min (a'), and so 
U—{U<) = exp(i<»(a')). 

Of course, due to gauge invariance, the fine configuration a M (x) solving 
Eq. ( |40|) is not unique; indeed, an entire orbit of solutions exists, all related by 
gauge transformations not altering the fixed coarse configuration a'^(x b ). It is 
possible to exploit this restricted freedom and fix the gauge in order to simplify 
the problem. 

We (partially) fix the gauge by requiring: 



a™ a (2x b ) = a™ n {2x b + fi) , /i=l,2. (A.l) 
In this gauge, the constraint: 

[a Al (2xfc) + a Al (2x i , + /t)] (mod 2tt) = a'^x b ) (A.2) 

has the solution: 



a™ n (2x b ) = a™ n (2x b + ft) = -ot'Jxb) + n^(x b )n , (A. 3) 

where n^(x b ) can be or ±1 (the sign depending on the sign of a^(x(,)); ex- 
ploiting the gauge freedom we can fix: n^{x b ) = 0. 

For a given lattice site of the blocked lattice x b , we define the block of the 
sites of the fine lattice belonging to x b as: {x : x — 2x b + Ail + A22, = 
0, 1, 2}. We see that in our particular gauge, the condition (|A.2[) fixes the eight 



fine links lying on the border of the above defined blocks (Eq. (A.3)). As a 
consequence, the residual minimizing problem decouples in the various blocks; 
indeed the remaining degrees of freedom are the links internal to the blocks 
(four for each block) , and links internal to different blocks do not communicate 
between themselves. 

Minimization inside each block can be now trivially performed. Here gauge 
invariance can be again exploited by fixing: 
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a'™(2x b + 2) = a™ a (2x b + 1 + 2) . (A.4) 

The result is: 



a™" 1 (2x6 + i) = -(c/ 2 (xb) + a' 3 (xb + i)) + (n^Xb) - n 2 (xi,)) w 
a'™{2x b + 2) = a™ a (2x b + 1 + 2) = 

1 - 7T 

-(a' (x b ) + a[(x b + 2)) + {n 3 (x b ) + n±(x b ) - n^Xb) - n 2 (x b )) — 
a 2 in (2x b + 1 + 2) = -(a' 2 (xb)+a' 2 (x b + i)) + (n^Xb) - n 3 (x b )) n , (A.5) 



where the integer numbers rii(x b ),i = 1,2,3,4, are involved in the consistency 
conditions 



F%(2x b )-n 1 (x b )2ir e (-7r,7r] , (A.6) 

and analogous relations respectively for F^ 2 (2xb + 1), F^ 2 (2x b + 1 + 2) and 
F^ 2 (2x b + 2), where F^ 2 (x) is the lattice held tensor for the solving fine con- 



figuration given by Eqs. ([A. 3D and (A.5). For this solution, one has 



F^(2x b ) -71,271 = F^ 2 {2x b + l) - n 2 2ix = F^(2x b + l + 2) - n 3 2n = 

1 7T 

F^(2x b + 2) - n 4 2ir = -F^ 2 '(x b ) - (n, + n 2 + n 3 + n 4 ) - , (A.7) 



where we indicate with F^ 2 '(x b ) the field tensor of the coarse configuration 
a'(x b ). Using the above relations, its is clear that the total action of the fine 
configuration, according to Manton's definition, reads: 



(mod27r)] 2 = 

X 

^J2[Fl' 2 '{x b )-{n 1 (x b ) + n 2 {x b )+n 3 (x b ) + n i (x b ))2iT} 2 . (A.8) 
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The absolute minimum of the r.h.s. of the above equation among variations of 
the m(xb)'s is given by: 



S min = ra*6)(mod27r)] 2 . (A.9) 



Observe that the set of integer numbers realizing the absolute minimum (A.9) 
satisfy in particular the consistency relation ( |A.6| ) and the other three analogous 
relations for F^{2x b + V), F^{2x b + 1 + 2) and F^{2x b + 2). Eq. ( |A~9| ) is the 
Manton action for the coarse configuration, apart from a factor 1/4, which 
renormalizes the coupling in a trivial way, (3 — ► f3/4. 

For the fine minimizing configuration the following relation holds (see Eq. 



F%(x) (mod27r) = - [F%(x b ) (mod27r)] , x = 2x b + AJ + A 2 2, A M = 0, 1 , 

(A.10) 

from which it is evident that the lattice topological charge operator (|47]) repro- 
duces itself under the RG transformation. 
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